Tower Top Daily Average
library(tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5 v purrr 0.3.4
## v tibble 3.1.5 v dplyr 1.0.7
## v tidyr 1.1.4 v stringr 1.4.0
## v readr 2.1.0 v forcats 0.5.1
## Warning: package 'readr' was built under R version 4.1.2
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag() masks stats::lag()
top <- read.csv('Tower_Top_0066_5306_0586_Full.csv')
top$Temp <- ((top$Temp - 32)*5)/9
top <- separate(top, Time, into = c('hours','minutes'))
top$hours <- as.numeric(top$hours)
top <- na.omit(top)
for (i in 1:nrow(top)) {
if (top$minutes[i] == 30) {
top$hours[i] <- top$hours[i] + 0.5
}
}
top2 <- top %>% group_by(Date, hours)
top2$Place <- "Top2"
boxplot(Temp ~ hours, data = top2, col = "salmon", main= "Boxplots of Tower Top Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Tower Middle Daily Average
middle <- read.csv('Tower_Middle_7433.csv')
middle <- separate(middle, Time, into = c('hours','minutes'))
middle$hours <- as.numeric(middle$hours)
middle$Temp <- ((middle$Temp - 32)*5)/9
for (i in 1:nrow(middle)) {
if (middle$minutes[i] == 30) {
middle$hours[i] <- middle$hours[i] + 0.5
}
}
middle2 <- middle %>% group_by(Date, hours)
middle2$Place <- "Middle2"
boxplot(Temp ~ hours, data = middle2, col = "green", main = "Boxplots of Tower Middle Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Tower Top + Tower Middle Comparison
library(plotly)
##
## Attaching package: 'plotly'
## The following object is masked from 'package:ggplot2':
##
## last_plot
## The following object is masked from 'package:stats':
##
## filter
## The following object is masked from 'package:graphics':
##
## layout
mis <- middle %>% arrange(Date, hours)
mis$Place <- "Middle"
for (i in 1:nrow(mis)) {
if (mis$minutes[i] == 45) {
mis$hours[i] <- mis$hours[i] + 0.5
}
}
mis2 <- mis %>% distinct(Date, hours, .keep_all = T) %>% group_by(hours) %>% mutate(Average = mean(Temp, na.rm = T)) %>% slice(1)
mis3 <- top %>% arrange(Date, hours)
mis3$Place <- "Top"
for (i in 1:nrow(mis3)) {
if (mis3$minutes[i] == 45) {
mis3$hours[i] <- mis3$hours[i] + 0.5
}
}
mis4 <- mis3 %>% distinct(Date, hours, .keep_all = T) %>% group_by(hours) %>% mutate(Average = mean(Temp, na.rm = T)) %>% slice(1)
mis4$Place <- "Top"
topmiddle <- rbind(mis2, mis4)
topmiddle$hours <- as.numeric(topmiddle$hours)
plot1 <- ggplot(topmiddle, aes(hours, Average, color = Place)) + geom_smooth(aes(x = hours, y = Average), method = "lm",
formula = y ~ poly(x, 23), se = FALSE) + geom_errorbar(alpha = 0.3, aes(ymin = Average - sd(Average, na.rm = T), ymax = Average + sd(Average, na.rm = T), width = 0.2))
plot1

Tower Base Daily Average
base <- read.csv('Tower Base 4377.csv')
base$Temp <- ((base$Temp - 32)*5)/9
base <- separate(base, Time, into = c('hours','minutes'))
base$hours <- as.numeric(base$hours)
base <- na.omit(base)
for (i in 1:nrow(base)) {
if (base$minutes[i] == 30) {
base$hours[i] <- base$hours[i] + 0.5
}
}
base2 <- base %>% group_by(Date, hours)
base2$Place <- "Base2"
boxplot(Temp ~ hours, data = base2, col = "salmon", main= "Boxplots of Tower Base Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Plots of Tower
mis <- middle %>% arrange(Date, hours)
mis$Place <- "Middle"
mis$minutes <- as.numeric(mis$minutes)
mis$hours <- as.numeric(mis$hours)
for (i in 1:nrow(mis)) {
if (mis$minutes[i] == 45) {
mis$hours[i] <- mis$hours[i] + 0.5
}
}
mis2 <- mis %>% distinct(Date, hours, .keep_all = T) %>% group_by(hours) %>% mutate(Average = mean(Temp, na.rm = T)) %>% slice(1)
mis3 <- top %>% arrange(Date, hours)
mis3$Place <- "Top"
mis3$minutes <- as.numeric(mis3$minutes)
mis3$hours <- as.numeric(mis3$hours)
for (i in 1:nrow(mis3)) {
if (mis3$minutes[i] == 45) {
mis3$hours[i] <- mis3$hours[i] + 0.5
}
}
mis4 <- mis3 %>% distinct(Date, hours, .keep_all = T) %>% group_by(hours) %>% mutate(Average = mean(Temp, na.rm = T)) %>% slice(1)
mis4$Place <- "Top"
mis5 <- base %>% arrange(Date, hours)
mis5$Place <- "Base"
mis5$minutes <- as.numeric(mis5$minutes)
mis5$hours <- as.numeric(mis5$hours)
for (i in 1:nrow(mis5)) {
if (mis5$minutes[i] == 45) {
mis5$hours[i] <- mis5$hours[i] + 0.5
}
}
mis6 <- mis5 %>% distinct(Date, hours, .keep_all = T) %>% group_by(hours) %>% mutate(Average = mean(Temp, na.rm = T)) %>% slice(1)
mis6$Place <- "Base"
topmiddle <- rbind(mis2, mis6, mis4)
topmiddle$hours <- as.numeric(topmiddle$hours)
ggplot(topmiddle, aes(hours, Average, color = Place)) + geom_smooth(aes(x = hours, y = Average), method = "lm", formula = y ~ poly(x, 23), se = FALSE) + geom_errorbar(alpha = 0.2, aes(ymin = Average - sd(Average, na.rm = T), ymax = Average + sd(Average, na.rm = T), width = 0.2)) + geom_rect(aes(xmin = 5.75, xmax = 6.25, ymin = 20, ymax = 30),
fill = "transparent", color = "black", size = 0.5, linetype = 2) + ggtitle("Plots of Tower")

ggplotly(plot1)
library(tidyverse)
MP150 <- read.csv("MP150_5354_7432.csv")
Lago0066 <- read.csv("Lago 0066.csv")
Lago5353 <- read.csv("Lago 5353.csv")
MP150 Plot Daily Average
MP150 <- read.csv("MP150_5354_7432.csv")
MP150 <- separate(MP150, Date...Time, into = c("dates", "Time"), sep = " ")
MP150$dates <- as.Date(MP150$dates, "%m/%d/%Y")
MP150 <- separate(MP150, Time, into = c('hours','minutes'))
MP150$hours <- as.numeric(MP150$hours)
MP1502 <- MP150 %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = MP1502, col = "salmon", main= "Boxplots of MP150 Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Lago Plot Daily Average
Lago0066 <- read.csv("Lago 0066.csv")
Lago5353 <- read.csv("Lago 5353.csv")
Lago5353 <- Lago5353[,-1]
colnames(Lago5353) <- c("Number", "Date...Time", "Temperature", "Humidity")
colnames(Lago0066) <- c("Number", "Date...Time", "Temperature", "Humidity")
Lago <- rbind(Lago0066, Lago5353)
Lago$Temperature <- ((Lago$Temperature - 32)*5)/9
Lago <- separate(Lago, Date...Time, into = c("dates", "Time"), sep = " ")
Lago$dates <- as.Date(Lago$dates, "%m/%d/%Y")
Lago <- separate(Lago, Time, into = c('hours','minutes'))
Lago$hours <- as.numeric(Lago$hours)
Lago2 <- Lago %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = Lago2, col = "salmon", main= "Boxplots of Lago Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

hello <- Lago2 %>% filter(dates == "2019-08-04")
Parahuaco Plot Daily Average
Parahuaco <- read.csv("Parahuaco RH.csv")
Parahuaco <- separate(Parahuaco, Date...Time, into = c("dates", "Time"), sep = " ")
Parahuaco$dates <- as.Date(Parahuaco$dates, "%m/%d/%Y")
Parahuaco <- separate(Parahuaco, Time, into = c('hours','minutes'))
Parahuaco$hours <- as.numeric(Parahuaco$hours)
boxplot(Temperature ~ hours, data = Parahuaco, col = "salmon", main= "Boxplots of Parahuaco Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

GP Plot Daily Average
GP <- read.csv("GP RH.csv")
GP <- GP[,-1]
GP <- separate(GP, Date...Time, into = c("dates", "Time"), sep = " ")
GP$dates <- as.Date(GP$dates, "%m/%d/%Y")
GP <- separate(GP, Time, into = c('hours','minutes'))
GP$hours <- as.numeric(GP$hours)
GP2 <- GP %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = GP, col = "salmon", main= "Boxplots of GP Temperatures", xlab = "Hours from Midnight (2019-07-26)", ylab = "Temperature in Celcius")

Murcialago (Rest) Daily Average
Murcialago <- read.csv("Murcialago5306.csv")
Murcialago <- separate(Murcialago, Date...Time, into = c("dates", "Time"), sep = " ")
Murcialago$dates <- as.Date(Murcialago$dates, "%m/%d/%Y")
Murcialago <- separate(Murcialago, Time, into = c('hours','minutes'))
Murcialago$hours <- as.numeric(Murcialago$hours)
Murcialago2 <- Murcialago %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = Murcialago2, col = "salmon", main= "Boxplots of Murcialago Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Parahuaco (Rest) Daily Average
Parahuaco <- read.csv("Parahuaco RH.csv")
Parahuaco <- separate(Parahuaco, Date...Time, into = c("dates", "Time"), sep = " ")
Parahuaco$dates <- as.Date(Parahuaco$dates, "%m/%d/%Y")
Parahuaco <- separate(Parahuaco, Time, into = c('hours','minutes'))
Parahuaco$hours <- as.numeric(Parahuaco$hours)
Parahuaco2 <- Parahuaco %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = Parahuaco2, col = "salmon", main= "Boxplots of Parahuaco Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Danta (Rest) Daily Average
dataset <- read.csv("Denta Region RH.csv")
dataset <- dataset[,-1]
colnames(dataset) <- c("Number", "Date...Time", "Temperature", "Humidity")
Danta <- separate(dataset, Date...Time, into = c("dates", "Time"), sep = " ")
Danta$dates <- as.Date(Danta$dates, "%m/%d/%Y")
Danta <- separate(Danta, Time, into = c('hours','minutes'))
Danta$hours <- as.numeric(Danta$hours)
Danta2 <- Danta %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = Danta2, col = "salmon", main= "Boxplots of Danta Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Titi Cabin (Sleep) Daily Average
dataset <- read.csv("TitiCabin9223.csv")
dataset <- separate(dataset, Date...Time, into = c("dates", "Time"), sep = " ")
dataset$dates <- as.Date(dataset$dates, "%m/%d/%Y")
dataset <- separate(dataset, Time, into = c('hours','minutes'))
dataset$hours <- as.numeric(dataset$hours)
dataset2 <- dataset %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = dataset2, col = "salmon", main= "Boxplots of Titi Cabin Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")

Titi Trail (Sleep) Daily Average
dataset <- read.csv("TitiTrail2115.csv")
dataset <- separate(dataset, Date...Time, into = c("dates", "Time"), sep = " ")
dataset$dates <- as.Date(dataset$dates, "%m/%d/%Y")
dataset <- separate(dataset, Time, into = c('hours','minutes'))
dataset$hours <- as.numeric(dataset$hours)
dataset2 <- dataset %>% group_by(dates, hours)
boxplot(Temperature ~ hours, data = dataset, col = "salmon", main= "Boxplots of Titi Trail Temperatures", xlab = "Hours from Midnight", ylab = "Temperature in Celcius")
